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The running of the Top quark mass is considered in the nonperturbative framework of the 
Schwinger-Dyson equation. Based on the input of physical pole mass meassured at the Tevatron the 
method provides the resulting mass function which is almost constant at low spacelike and timelike 
scales. The skeleton loops including Standard Model Higgs and gluons are taken into account. The 
dominant two-loop skeleton contribution with triplet Higgs interaction is considered in addition to 
one loop dressed approximation of the top quark self-energy. 
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I. INTRODUCTION 



The quark masses are fundamental parameters of the Standard Model. The precise knowledge of quark masses at 
various scales is important for several reasons. In hadronic physics it is necessary for precise determination of CKM 
matrix elements, while the theoretically extracted information about values of quark masses at very high momenta 
can be useful for model builders. 

In perturbation theory QCD approach the definition of the running mass is based on the renormgroup evolution 
equations (RGEs). MS bar scheme represents short range mass definition and is commonly used due to its technical 
simplicity [l], Q ■ On the other side the RGEs method cannot provide reliable results at low momenta where the 
perturbation method fails. A legitimate question is what is the relevant scale of the applicability of the perturbation 
theory, when the corrections to the quark masses are evaluated. To determine this, the running masses calculation 
has to rely on nonperturbative QCD techniques. So far, there are two methods that directly follow from the first 
principles: the first is lattice theory (for a review see |, Hi, 1,0,1) which is based on the discretized Euclidean 
space, the second is the functional method represented by a continuous framework of QCD Schwinger-Dyson equations 
[{J, [l(| • Within certain phenomenological assumptions the QCD sum rules [ill, G3 are used to determine the quark 
masses. 

After the top quark discovery, the top quark mass value is obtained in the fairly limited regime, the CDF fl3j 
and DO [T3| collaborations measure the resonant top quark mass. The particle data group [l5[ quoted the value 
M t = 174.2 ± 3.3GeV as the pole mass of the top quark. 

The evolution of the Yukawa coupling has already been studied in [l[ . However, the Yukawa interaction has not been 
considered in RGEs for the top quark running mass. The contribution to the quark selfenergy due to the Higgs boson 
has been studied in SDE framework for the first time in (lfjj . This study has been performed with two approximately 
equivalent inputs: M(0) = 179GeV and M(Mf) = 174GeV, noting that the later was defined incorrectly at spacelike 
scale M t . In the present paper we go beyond the one loop approximation and take into account the two loop Higgs 
contribution as well. 

There is a striking evidence that the RGE perturbation calculation overestimated low scale top quark mass from 
the very beginning (going from high spacelike Q 2 to the infrared values) . Recall that in the paper 1] the renormgroup 
equation for top quark running mass M{p 2 ) has been solved in MS bar renormalization scheme with the following 
result (in GeV): 

M(180 2 ) = 170.1; M (M 2 ) = 170.8; M(91.2 2 ) = 180; 

Af(4.3 2 ) = 253; M(1.3 2 ) = 318; M(l) = 339 , (1.1) 

for spacelike arguments in the brackets and have been quoted within ~ (12 — 25)GeV error due to experimentally 
determined physical mass ( to that date, it was M t — 180GeV). 

Recall that the physical pole mass M t is determined in Minkowski space as the S (Mt) = 0, in other words 
M t (—Mf) = M t . Assuming that the fit procedure of M(M 2 ) from M t used in [l[ and developed originally in [TtJ is 
reliable (note, the relation between MS mass and on shell mass is recently known to the order a 3 ), one necessarily 
must conclude that the renormgroup evaluation of masses becomes unreliably overestimated below the scale /1 ~ M z . 
While for leptons and light quarks perturbation QCD works perfectly at the M z scale, it appears that for an accurate 
estimate of the running top quark mass at M z mass scale might not be adequate. Technically this is because already 
one loop correction to M is enhanced like 



5M ~ olqcdM . 



(1.2) 



FIG. 1: Diagrammatic representation of top quark SDE. The solid line stands for quark and the dot (dashed) line stands for 
the gluon (Higgs) respectively. The circles represent full vertices. 



In other words, the exceptionally large mass of the top quark itself spoils the usual correctness of perturbative QCD 
at electroweak scale. 

This is one of the main reason of the present study to calculate the evolution of top quark mass in the whole 
momentum of range, obtaining thus correct information for the low energy scales. In perturbative MS schemes the 
running mass grows from MS value m(m) = 170 to fi — Mz about amount of 10GEV and further blows up when 
evolved to the infrared. We will argue that M(M t ) and M{Mz) differs about tiny amount ~ l — 2GeV and the running 
top quark mass function remains stable when using selfconsistent framework of our SDE equations. The knowledge 
of here observed infrared stability of the top quark mass should be useful whenever a selfconsistent treatment is 
required, i.e. for instance when one considers Higgsonia fl8l . Il9j and the effect of top quark loop in the equations for 
Higgsonium bound states. Further, the top quark circulates in the loop of penguin diagrams describing rare mesonic 
electroweak decays (see e.g. 0, Hl|)- In this case the typical energy of decaying B mesons is of the order nib, so the 
good knowledge of the infrared value of the top quark mass is important for the description of heavy meson decays. Of 
course, the knowledge of mass at high scales can useful for model builders. However, in the nonperturbative treatment 
here we are mainly for the physic not far above the electroweak scale, the knowledge of the quark mass at higher 
scales can be useful for model builders as well. 

The last but not least motivation is a direct check the effect of higher order corrections including Higgs trilinear 
coupling on the solution. To do this the appropriate two loop skeleton diagram is calculated and included into the 
top quark SDE. These, and other details of the model are described in the Section II. and Section III. of presented 
paper. To find the correct solution in the full Minkowski space is a problematic task for a strong coupling theory like 
QCD. First we solve quark SDE in Euclidean space by standard numerical manner in the Section IV. In Section V. 
we continue the solution to the timelike axis in a way that experimentally known pole mass is achieved by the correct 
solution. This is achieved by resolving of the SDE with renormalized mass adjusted to obtain the correct physical 
pole mass at the end. 



II. SCHWINGER-DYSON EQUATION FOR TOP QUARK MASS FUNCTION 

Neglecting the weak interaction, the quark propagator S can be conventionally characterized by two independent 
scalars, the mass function M and the renormalization wave function Z such that 

The SDE for the inverse of S reads 

Sip)' 1 = i>-9Y <^>-E A (p)-E h (p) + ... (2.2) 

Sfc(p) = igy J jr^4 T h{q,p)G h (p - q)S(q) 

where gy is the top Yukawa coupling, Higgs vev < <j> >= 246GeV^/\/2 and g is QCD gauge coupling. The dots 
represent omitted contributions, e.g. W, Z, 7 and related Goldstone exchanges. G, V stand for boson propagators and 
the top quark-boson vertices and they satisfy their own SDEs. 

The knowledge about these Greens function is necessarily limited due to theoretical and experimental reasons. They 
need to be approximated if they are not selfconsistently contained in a given truncation scheme of the SDEs system. 
A natural treatment of this problem is to make an expansion in the number of loops. Performing such an expansions 
for vertices T = '}2 li T l and substituting this into the selfenergy ()2.2|) . one gets the expansion for the mass function. 
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Explicitly, the loop expansion for the selfenergy in (|2.2[) should read 

^ a = ^ ] +J2^a +1] (2-3) 

i 

and similarly for £/,. 

In the simplest approximation the first order estimates can be obtained by using the classical vertices 

E£(p) - i9v J ^G h (p-q)S(q) (2.4) 

where the all propagator functions entering the Eqs. (|2.4[) are fully dressed. 

Including "radiative corrections" to the SM model Higgs one should get coupled SDEs for G and S. In the case of 
light Higgs, the top-antitop quark loop contribution would lead to the extremely large negative contribution to the 
Higgs boson mass. This mass hierarchy problem, although formally solved by renormalization, is one of the main 
motivation for extension of the Standard Model and the reason why the SM is regarded as an effective low energy 
theory. In the extensions of SM the mass hierarchy is stabilized by the introduction of the other scalars [2J, [25j, [26j , 
SM doublets [UHl], or is eliminated by supersymmetry or the Lee- Wick SM modification (29|. In all these models, 
a new particle content is expected at few TeV, the quadratic divergences to Higgs mass are reduced and the free 
propagator could be a reasonable approximation of the exact Higgs propagator for a broad regime of scales. Therefore 
the simplest -free Higgs boson propagator: 

G h (p) = — L- ^ (2.5) 



is used, where m is the physical Higgs boson mass (|2.6|) . 

Following the recent precision test of the Standard Model [22} ■ the analysis of the radiative corrections favor a light 
Higgs boson m ~ 76 GeV. Because of the lack of an experimentally observed Higgs particle, the mass of the Higgs 
boson could be rather close to the experimental lower bound m > 114.4 GeV (23|. In this paper the following value 
of the Higgs mass is chosen 

m = 120 GeV, (2.6) 

as the input parameter in our model. 

At low scales, q ~ Aqcd, the running QCD coupling is large and the dressing of the gluon-quark-antiquark vertex 
can play an important role in the description of light flavor dynamics [30j ]. However, in the case of the top quark, the 
running coupling becomes quite small cxqcd(M) —0.1 and one can economically include the contribution of higher 
orders to the effective running coupling. For this purpose the following prescription for the SDE QCD-part kernel is 
used: 

g 2 G^(k)T u (q,p) -> 47ra(fc 2 , A) ~ g fc2 + gH ^ . ( 2 .7) 

where a represents the analytical running coupling [32], l33l . l34l |35| . In the one loop approximation it is given by 
the following expression: 

a(q Aqcd) = / du; 2 — , (2.8) 



where 



o q - v 



4tt//3 



Pg {u, A Q cd) = 2 TT /A2 r ■ ( 2 - 9 ) 

7r 2 -ln (uj/A qcd ) 

Recall that the analytical running coupling is constructed in a simple way that avoids the unwanted artifact of 
perturbation theory- the Landau pole at q 2 = Aq CD - which is subtracted away and thus the running coupling is 
free of unphysical singularities. The procedure has been generalized to higher orders, provided that the ultraviolet 
asymptotic behaviour of such running constant is identical with the perturbative result. In the numeric here the one 
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FIG. 2: Diagrammatic representation of selfenergy contribution due to the Higgs. The circle stands for the full vertex, the 
rhs. is the approximation employed here. 

loop approximation (|2.9p is used with the numerical value of Aqcd, ^qcd — 500MeV for six active quarks. The 
beta function coefficient is thus 

UN C -2N f , 
P= ^ L (2.10) 

with iV c = 3, N f = 6. 

The computation is carried out in Landau gauge and the Z = 1 approximation is used. Whilst in pure gauge theory 
the effect of the Z = 1 approximation can be minimized by proper adjustment of the gauge fixing, the importance of 
Z in the presence of the Higgs field is not explored and remains to be estimated in a future study. 

III. SOLVING TOP QUARK SDE IN EUCLIDEAN SPACE 

Using the following formula 

dr\ sin 2r q 



7r I 2 — 2\l\ \p\cosrj + p 2 + m 2 



_p2 _ ^2 _ m 2 _|_ ^/(p 2 + q 2 + to 2 ) 2 — 4p 2 q 2 

' ( } 

the angular integrations in one loop skeleton diagram in (|2.4p can be easily evaluated. After the explicit integration 
the Higgs-top loop contribution can be cast into the one dimensional integral 

4 1] (*) = ^ rdy-^K( X ,y, m >), (3.2) 
7T J y + M 2 {y) 

where ay — gy/^ 71 an d the functions K is defined as 

K(x,y,z) = ^ = . (3-3) 

x + y + z+ ^(x + y + z) 2 — Axy 

Likewise, for the one loop QCD contribution we get 

S W (x) = - dy ™M V(x, y) , (3.4) 
irj y + M 2 {y) 

V(x, y) = - [°° duo^- [K(x, y, 0) - K{x, y, to)} . (3.5) 
Jo u 

In addition, the one loop skeleton contribution to the Higgs-quark-antiquark proper vertex (see Fig. [2]) is carefully 
included. This is equivalent to the two loop 1PI contribution for the top quark dynamical mass function which reads: 

4V) = A™#I(p), (3.6) 
where A/4 is the Standard Model quartic Higgs coupling and I(p) is the following two loop integral: 



with the function V defined as 



. . Tr f dH 1 I+MO) 

Hp) = -r* 



X I 



4 J (2tt) 4 (I - p) 2 - m 2 I 2 - M(l) 2 
d 4 q 1 A + M(q) 



(2tt) 4 (q - p) 2 - m 2 q 2 - M{q) 2 
(q - I) 2 - m 2 ' (3 ' 7 ^ 
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After the Wick rotation to Euclidean space, the integrations in (|3.7p are not calculable directly, however most of them 
can be calculated analytically by performing just one quite standard angular approximation (|3.7p ). This approximation 
eliminates the angle between the two internal loop momenta in the following manner: 

(I - qf -> l 2 6(l 2 - q 2 ) + q 2 6{q 2 - I 2 ) (3.8) 

and so writing also for l.q product (this stem from Dirac trace) 

l 2 + q 2 -{l-qf 



l-q = 



q2 8(l 2 -q 2 ) + l le(q 2 -l 2 ), (3.9) 



the expression for / can be recast as: 

m = 



2 v 1 ' 2 
dH 1 f d 4 q 1 



(2tt) 4 [I - p) 2 + m 2 J (2tt) 4 (q - p) 2 + m 2 
M {l) M i+q) y,2 e{l 2 q2) + M}MM+tll e{q2 p) 



{q 2 +M 2 (q))(l 2 + M 2 (l)) 

Here, it is an opportune point to remark that such an angular approximation has been extensively used in phenomeno- 
logical SDE studies of QCD and QED4 even at one l oop level. In our case the coupling constant is small enough and 
following the critical one loop analysis performed in [3l| , this must be a reliable approximation in our two loop case. 
In Euclidean domain it can lead to a few percent error in I. As we have estimated posterior, it entails only a tiny (a 
few promile) error in the total result for M. 

Using the formula (|3.1[) the remaining angular integrations can be performed, resulting the following expression for 

I: 



, . dq 2 p 2 + q 2 + m 2 — \/(p 2 + q 2 + m 2 ) 2 - 4p 2 q 2 

W = ' 8^ P 2 [g 2 + Af 2 (g)] 



dl 2 p 2 + l 2 +m 2 - ^(p 2 + I 2 + m 2 ) 2 - 4p 2 l 2 
8^ p 2 [l 2 + M 2 {l)] 

M (l) MMW /2 e{l2 q2) + M (l) M^,2 9{q 2 _ g 
X (q 2 +M 2 (q))(l 2 + M(l) 2 ) 

(3.10) 

In what follows we interchange of the variables / <-> q in the second term of the third line of the Eq. (|3.10p . Considering 
the appropriate prefactors, I(p) can be finally written in the following way: 



i , 2 P 2 + q 2 + m 2 - J(p 2 + q 2 + m 2 ) 2 - 4p 2 q 2 

m = L dq s^i 



f 9 2 p 2 + l 2 +m 2 - y/(p 2 + l 2 + m 2 ) 2 - ApW 
; ' p 2 [l 2 + M 2 (l)] 

_L / 2 /2 + M(l)M(g) 

X 32tt 2 (q 2 + m 2 ){q 2 + M 2 (q))(l 2 + M 2 {1)) 



(3.11) 



Putting these all together, the SDE for top quark mass function that is to be solved reads 

M(x) =g Y < > +S^ ] (x) + ^)}{x) + Y, l2] (x), (3.12) 

where the individual terms are given by (|3.4|) . (|3.2p and (|3.6p wherein / is given by Rel. (|3. 1 1|) . As a consequence of 
the Z = 1 approximation the function M{x) is renormgroup invariant. After making a subtraction the "renormalizcd" 
equation actually solved reads 

M(x) = M(0 + £(*) - E(0 

S(.x) = sW(x) + sW(x) + sl 2] (x), (3.13) 

where the renormalizcd mass at the scale ^ is related to the bare top quark mass through the following rel.: M(£) = 
g Y < (j)> +£(£). 
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FIG. 3: Running Top quark mass as described in the text. The solid line represents the full solution, dot-dot-dashed line 
stands for the case when two loop skeleton is omitted. The dashed and dotted lines stand for pure QCD and Yukawa solutions 
respectively. 

IV. RESULTS IN SPACELIKE REGIME 

In this section, we discuss the numerical solution of the renormalizcd Euclidean SDE (|3.13[) . The physical mass 
pole being on the timelike axis cannot be directly used for the solution. The main purpose of this section is to exhibit 
the importance of various contribution in the case of light Higgs exchanges. 

The SDE (|3.13[) has been solved by the method of iterations with high accuracy. For this purpose we have chosen 
the (spacelike) renormalization scale to be 

f = IQQGeV 2 (4.1) 

and fixed the renormalizcd mass M(£) through the Yukawa coupling. 

The two loop diagram depicted in Fig. [^includes the triplet Higgs interaction constant, which is already determined 
through the quartic one. In our numerical calculation the coupling constant actually used is read from the relation 
A = (at the given scale £). 

The resulting mass function is displayed in the Fig. [3] The presented calculations were performed with gy — 0.977, 
so the corresponding renormalized mass is adjusted so that M(100GeV) — 169.25GeV. With these inputs we get the 
numerical solution. As the presented solution is regularization independent, we used hard cutoff regulator A >> M(0) 
obtaining the same solution when A was varied through many orders. The mass function is increasing when going to 
infrared, reaching its infrared value A/(0) = \7Q.QQGeV , being not far the experimental one. How to gain the solution 
actually based on determined physical top quark mass will be discussed in the section. Before this we discus some 
general features of the solution. 

In our presented framework of SDEs the dynamical mass function is slowly varying function in the infrared. Up to 
few GeV contribution the infrared mass does not change drastically at the scale of 0-100 GeV. 

The Yukawa interaction between Higgs and top quark is quite strong even when comparing to the QCD interaction 
strength. In Fig. [3l we show the comparison of solutions stemming purely from the Yukawa interaction and from 
the pure QCD (by switching off QCD or Yukawa interaction). The same value of the renormalized top quark mass is 
kept for this purpose. As expected, the QCD dominates in the infrared regime, while in high momenta, q 2 >> M 2 
both interactions are of the same magnitude. 

Interestingly, the two loop effect with Higgs trilinear coupling gives a marginal contribution for all p 2 . Numerically, 
two loop skeleton effect is comparable with the one loop PT electroweak corrections. For a heavier Higgs the one loop 
Higgs contribution becomes less important, while the two loop contribution appears to be less affected since the triplet 
Higgs coupling is getting strong. We have also solved the SDE with different Higgs masses as well. For instance Higgs 
heavy as m — O.bTeV, two loop contribution becomes more important giving rise a few GeV negative contribution in 
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the infrared top quark mass. However, one should note that in this case the Higgs sector becomes strongly interacting 
what would require more careful reinvestigation due to the new nonperturbative dynamics [42l l43l |44| . 

V. SOLUTION FOR ALL MOMENTA 

Experimentally the top quark mass is reconstructed by collecting jets and leptons. From the position of the bunch 
in cross section measured at the Tevatron the pole position is identified M t — f72.4 ± 1.4. The ambiguity and 
uncertainty of the full top propagator pole mass is affected by experimental methods and theoretical weaknesses of 
perturbation theory description of jets, e.g. reconciling the contribution of soft and collinear particles. Furthermore, 
the correct identification of the mass requires nonperturbative technique at all. While, including perturbative 1-loop 
6, W electroweak correction this pole could only move into the second sheet complex plane giving rise top quark decay 
width T t = l.hGeV, the perturbation theory cannot give nonambiguous result because of uncertainty proportional 
to Aqcd [H, H3, [Hj]. The real pole of the (pure QCD) perturbation theory can turn to be complex because of 
confinement phenomena as recently observed in [39| by studying complex mass generation in temporal Euclidean 
space. 

In this paper we do not solve the problem of confinement of th top quark in SDE framework, instead we show 
that the running top quark mass function turns to be stable, very slightly varying, quantity when continued to the 
timelike momenta. For this purpose the mass function M{— x) at timelike square of the fourmomenta t — —x > 
is constructed by the continuation of the left hand side of SDE by taking x — > —t. We can write for the continued 
solution 

M{t) = M(0 + £(«)-£(£) (5.1) 
£(t) - £™(t) + £W(t) + £L 21 (*). 

^^Vf^^,^), (5-2) 



K(-t,y,z) = 



'2y 



-t + y + z+ ^(-t + y + z) 2 + 4ty 



and similarly, the functions £^ , £^ are obtained by the substitution x — > — t in their kernels. 



Since the mass function on the rhs. of Eq. (|5.2|) remains defined at the spacelike regime, the pole mass cannot 
be used as the renormalized point directly. To achieve the solution of SDE with experimentally known value of top 
quark mass, we shift the renormalized mass in (|3 . 1 3|) and then have a look for the solution for M t by integrating the 
equation (|5.ip . With sufficient accuracy it is easily achieved by hand iteration process. 

The experimentally observed mass knowledge based solution is presented in Fig. [U The numerical value M t = 
172.6GeV is obtained as the solution for pole mass. The timelike solution is plotted at the negative axis. The resulting 
Yukawa coupling to our 12QGeV heavy Higgs field has been adjusted as <?y = 0.9845. The solution is real everywhere 
as the method is inefficient to provide absorptive part from the nonanalytical cut at real axis at M t . The experimental 
uncertainty defines the errors repesented by narrow band of width ~ 2 — 2.8GeV with presented solution inside. We 
do not display these. 

The other interesting values we can quote here are (in GeV): 

M(10 12 ) = 134.5; M(M t 2 ) = 171.0; M(M|) = 172.4; M(0) = 173.0; M{-M 2 ) = 172.6, M(-10 12 ) = 151.2, . (5.3) 



VI. CONCLUSION 

The SDE calculation of running top quark mass previously discussed in the literature [l6[ is presented in some 
extent. The obtained solution is based on the measured physical top quark mass. The top quark mass can be safely 
evolved to small q 2 when one avoids the pathology of perturbation theory, e.g. Landau pole in gluon propagator. It 
exhibit great stability at all scales of spacelike and timelike domain as well. For the timelike domain the function is 
such slowly varied that the top quark physical mass appears to be a rather good approximation at all low scales. 

At low scales, QCD contribution dominates over the one due to the Higgs loop(s), at large q 2 both Higgs and QCD 
loops are comparable. In addition, we estimated the two loop Higgs contribution, which gives only tiny contribution 
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FIG. 4: Running top quark mass based on observed Tevatron pole mass. The solid line represents the spacelike and dashed 
the timelike solution, dotted solid line represents the linear function / = \f[x) ,which when cuts the dashed line, identifies the 
real pole of the full top quark propagator). 



for the case of the light Higgs. The extension of presented technique to the more general models, e.g. with more Higgs 
doublets and/or scalar singlets added to SM Higgs sector is straightforward. 
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